The modular /S-matrix as order parameter for 
topological phase transitions 



F A Bais 1,2 and J C Romers 1 

1 Institute for Theoretical Physics, University of Amsterdam, Science Park 904, P.O.Box 94485, 
1090 GL Amsterdam, The Netherlands 

E-mail: j . c . romer s@uva . nl 

2 Santa Fe Institute, Santa Fe, NM 87501, USA 

Abstract. We study topological phase transitions in discrete gauge theories in two spatial 
dimensions induced by the formation of a Bose condensate. We analyse a general class of euclidean 
lattice actions for these theories which contain one coupling constant for each conjugacy class of 
the gauge group. To probe the phase structure we use a complete set of open and closed anyonic 
string operators. The open strings allow one to determine the particle content of the condensate, 
whereas the closed strings enable us to determine the matrix elements of the modular S-matrix, 
both in the unbroken and broken phases. From the measured broken S'-matrix we may read off 
the sectors that split or get identified in the broken phase, as well as the sectors that are confined. 
In this sense the modular S'-matrix can be employed as a matrix valued non-local order parameter 
from which the low-energy effective theories that occur in different regions of parameter space can 
be fully determined. 

To verify our predictions we studied a non-abelian anyon model based on the quaternion group 
H = D2 of order eight by Monte Carlo simulation. We probe part of the phase diagram for the 
pure gauge theory and find a variety of phases with magnetic condensates leading to various 
forms of (partial) confinement in complete agreement with the algebraic breaking analysis. Also 
the order of various transitions is established. 



PACS numbers: 73.43.-f, 71.10.Pm 
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1. Topological order and topological symmetry breaking 

1.1. Introduction 

The study of hidden symmetries and phase structure of gauge theories has a rich 
history A central motivation in the early work was the understanding the confinement 
phenomenon in non-abelian gauge theories. An important step forward were the lattice 
action formulation by Wilson pQ and its Hamiltonian version by Kogut and Susskind 
[2] , which allowed for an expansion around the strong coupling limit where confinement 
is manifest. The question was whether the confinement regime would extend all the 
way down to zero coupling. Wilson also introduced his celebrated loop operator as a 
diagnostic for confinement in the pure gauge theory. 

Not long thereafter 't Hooft and Mandelstam suggested confinement in 3+1 
dimensions to be a consequence of the dual Meissner effect, caused by magnetic disorder, 
notably a condensate of topological degrees of freedom, be it monopoles or fluxes [3J. 
't Hooft explained the (2 + l)-dimensional version of confinement in SU(N) gauge 
theories as a consequence of the condensation of Z N fluxes corresponding to the center 
elements of the gauge group [3J. Polyakov on the other hand proved confinement of 
(2 + l)-dimensional compact QED due to monopoles |5] and showed furthermore that 
the finite temperature deconfinement transistion in d = 3 + 1 is due to a Wilson 
line (string) condensate [6]. The interpretation of such transitions is related to the 
spectrum of topological defects in the various Higgs phases. For example for U(l) in 
3+1 dimensions the occurence of a phase transition depends the compactness of the 
group, i.e. on the presence of monopoles. For non-abelian groups the situation depends 
on the representation of the Higgs field. In the case of SU(2) for example one may or may 
not find a phase transition between the confined and Higgs phases depending on whether 
the breaking is achieved with two adjoints or a fundamental representation. Shenker and 
Fradkin [7] constructed the phase diagrams for the SU(2) (lattice) gauge theory for the 
distinct cases with Higgs fields in the adjoint and the fundamental representation, and 
proved the absence of a transition between confined and Higgs phase in the latter. The 
difference can indeed be traced back to the topological structure: in the broken phase 
with a fundamental Higgs no defects exist because the group is broken completely. In 
the broken phase with one isovector, there is a residual U(l) and one has monopoles, 
and these will be confined if one adds another generic isovector. However with the two 
vectors there is still the Z 2 center of the group which is unbroken, which implies the 
existence of Z 2 fluxes in the broken phase. 

Lattice formulations of the Z 2 theories go back till the early seventies [HI E] and were 
also studied in [7] , where also the phase diagram of the Z 2 gauge theory with a matter 
field was computed in detail. A Hamiltonian formulation of Zjy gauge theories and 
their duality properties in two and three dimensions was made in pi)], where the same 
algebraic structure ([T| was constructed out of generalized Ising type lattice models. All 
authors agree that the transition to the confined phase is due to the condensation of 
magnetic fluxes. The analysis is facilitated by the introduction of nonlocal order and 
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disorder parameters such as the Wilson and 't Hooft loops and the algebra they form, 
from which the phase structure could be understood qualitatively. As these operators are 
nonlocal, the algebraic structure exhibits the nontrivial braid proper ies that characterize 
the phase. For Zjv theories one obtains an underlying ZJ^ ® Z^ aff symmetry of which the 
electric, magnetic and dyonic (anyonic) charges form irreducible representations (n, m) 
labeled by a pair of integers mod iV (m, n = 0, 1, 2, . . . , N — 1) . If one defines the loop 
operators A^ n ' m \C), then 't Hooft derived the crucial algebraic relation [3J: 

^(n,m) ^( n 'i m ') (C') ^( n 'i m ') (C') A^ n ' m ^(C) g27ri(nm'+mn')fc/A'' /-|\ 

where k e Z is the linking number of the loops C and C . These loops can be interpreted 
as linked timelike loops of worldlines obtained after creating and subsequent annihilation 
of two particle-antiparticle pairs. But they can also be interpreted as two intersecting 
spacelike loops on a torus winding around the two different cycles, one being a Wilson 
line and the other a magnetic Dirac string. On the other hand one may think of equation 
as specifying the braiding relations for the representation theory of a Z^ x 7h N 
algebra endowed with a (unique) nontrivial braiding structure that exactly belongs to 
the quantum double algebra of Z^ denoted as D(Zn). This operation on the torus is 
a topologically nontrivial vacuum to vacuum transition (because no charges or fluxes 
are left) and that shows that there is a periodic vacuum structure in the theory, leading 
to a an iV 2 -fold vacuum degeneracy on the torus. This degeneracy in turn equals the 
number of particle species, i.e. the total number of unitary irreducible representations 
of -D(Zjv), which for the case at hand equals N 2 . Similar algebras have surfaced in 
the study of emergent gauge theories from closely related lattice models like Kitaev's 
toric code [H] and the Levin- Wen type spin models p2] , where furthermore the vacuum 
degeneracy (of a gapped phase) on the torus is proposed as a criterion for topological 
order. 

It was also early on pointed out that "the" Higgs phase structure of a gauge theory, 
where no massles gauge degrees of freedom (no continuous groups) survive may be very 
rich by itself, allowing for distinct phases which exhibited different spectra of purely 
topological degrees of freedom [T3]. By breaking with sufficiently high-dimensional 
representations of SU (2) for example, one finds that not just the center, but any discrete 
subgroup of SU(2) can be selected to survive, also non-abelian groups like D N , or the 
tetrahedral group etc. In the corresponding phases one finds a rich variety of non-abelian 
fluxes that exhibit highly nontrivial braiding and fusion properties with charges and also 
among them selves. The effective low energy theory in such a Higgs phase is also called 
a discrete gauge theory (DGT), these are indeed gapped phases with only topological 
degrees of freedom, and therefore by definition phases with topological order. The full 
(2 + l)-dimensional description of all the anyonic sectors of discrete gauge theories was 
given in [H] and revealed an underlying hidden quantum group structure, i.e. the 
excitations of a discrete gauge theory with gauge group H, were shown to correspond to 
the irreducible unitary representations of the quantum double D(H) of H as defined by 
by [IS]- The most interesting aspect of this perspective being that the ordinary "electric 
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charge" and topological "magnetic flux" degrees of freedom are both part of the same 
representation theory, the magnetic degrees of freedom labeled by the conjugacy classes 
of H and the electric parts by the representations of the centralizer group in H of 
the given magnetic flux. This leads to a rather intricate interdependence of electric 
and magnetic sectors in case the goup H is non-abelian and a different way to look at 
order and disorder operators. We return to D(H) and its representations in Section 



2.2 An important consequence of this underlying symmetry is that it suggested that 
in the pure DGT a corresponding set of gauge invariant loop operators should exist, 
thereby generalizing the Wilson/'t Hooft operators to all anyonic species, and indeed 
such operators have been constructed [IB]- With these operators in hand one can now 
probe the full phase structure of discrete gauge theories in a lattice formulation and 
that is what will be done in this paper. 

The hidden quantum group symmetry and its representation theory which forms 
a modular tensor category, is a powerful way to characterize the distinct topological 
phases of pure DGTs and forms the natural basis for the concept of topological symmetry 
breaking, which is defined as the breaking of quantum group symmetry [171 HH [19J. A 
brief review of this symmetry breaking mechanism is given in section 1.3 In this paper 
we study the validity of this mechanism by probing the various groundstates through 
measuring the expectation values of single as well as linked anyonic loop operators, as is 
explained in Section L4 In recent years we have witnessed a growing interest in systems 
that allow for the realisation of different topological phases, the examples can be found 
in Levin- Wen models [20], the Kitaev Honeycomb model [21], discrete gauge theories 
[17] . and last but not least in quantum Hall systems [221 E31 EH [251 EE]. Most of these 
can be understood from the point of view of topological symmetry breaking involving 
the formation of a Bose condensate. In most instances a Hamiltonian framework is used, 
but in this paper we show that to analyse such systems it may be profitable to switch 
to a euclidean action formulation which allows for the use of Monte Carlo simulations 
to determine the phase structure. In this approach it is easy to directly measure the 
modular ^-matrix in any phase of the system which explains why we call this matrix 
an orderparameter of such a system. After briefly recalling the basic ingredients of 
topological order and topological symmetry breaking and settle some notation, we show 
how open string operators and the modular S-matrix can be used as order parameters 
and phase indicators for topological symmetry breaking. We then introduce a class of 
multi-parameter lattice actions for non abelian discrete gauge theories and verify the 
theoretical analysis [181 US] in detail by numerical simulations. 



1.2. TQFT basics 

In this section we set the stage and fix the notation for the rest of this paper. We 
study phases of systems that are described by a Topological Quantum Field Theory 
(TQFT) in 2 + 1 dimensions. We label the different sectors or (anyonic) particle species 
by a, b, c, . . . . The two interactions between two particles in a TQFT are fusion and 
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braiding. 



Fusion We describe fusion by the rule 



a x b = N f c , 



(2) 



c 



where the integer multiplicities iV" give the number of times c appears in the fusion 
product of a and b. The fusion algebra is associative and commutative, and has a unique 
identity element denoted as "1" that represents the vacuum. Each sector a has a unique 
conjugate a (representing the corresponding anti anyon) with the property that their 
fusion product contains the identity: 



Braiding The particles in a 2 + 1 dimensional TQFT can have fractional spin and 
statistics. Rotating a particle a by 2tc (also called twisting) multiplies the state vector 
by a phase equal to the spin factor 9 a 

i \ twist p. I \ 

| a) ->■ 9 a \a), 

generalizing the usual +1 (—1) known from bosons (fermions) in 3 + 1 dimensions. 
Adiabatically moving a particle a around another particle b in a channel c is called a 
braiding and can have a nontrivial effect on the state vector of the system, given by 



Quantum dimensions The quantum dimensions d a of particle species a are another 
set of important quantities in a TQFT. These numbers satisfy the fusion rules (|2]), i.e. 
d a db = J2 c N^ b d c . The quantum dimension of an anyonic species is a measure for the 
effective number of degrees of freedom, corresponding to the internal Hilbert space of the 
corresponding particle type. The Hilbert space dimension of a system with N identical 
particles of type a grows as (d a ) N for N large. In general, the quantum dimensions d a 
will be real numbers; however for DGTs they are integers. The total quantum dimension 
T> of the theory is given by 



and the topological entanglement entropy of the ground state is proportional to log I?. 



a x a = 1 + 




9 c /9 a 9b- 




Diagrammatics There is a powerful diagrammatic language to express the equations 
describing the TQFT, which we will use to relate the values of observables as they can be 
measured in the different phases. In this paper we will use the notation and definitions 
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given by Bonderson [27]. Particle species are represented by lines, fusion and splitting 
by vertices. A twist is represented by a left or right twist on a particle line: 




The evaluation of simple diagrams is rather straightforward, and complicated 
diagrams can be simplified using braid relations and the socalled F symbols which 
follow from associativity of the fusion algebra. The simplest examples are the closed 
loop of type a that evaluates to the quantum dimension d a : 

°Q=d a , (4) 

whereas the twisted loop equals d a 9 a : 



a d a (5) 



Of particular interest are the generators of the modular group, S ao 

. s* = 3Z 1 *±i« (6) 

u. b c 

and T a b — e -27 ™^ 24 ^^,^ where the c is the central charge of the theory, not to be 
confused with a particle type. 

As we mentioned before the importance of the rather abstract diagrammatic 
notation is that the diagrams directly correspond to observables in our euclidean lattice 
gauge theory formulation. In the euclidean three dimensional formulation of topological 
theories the values these diagrams have, correspond to the vacuum expectation values 
of the corresponding anyon loop operators, for example in the unbroken phase one may 
measure 

0)o = (?) 
where the LHS is now defined as the value of the path integral with the nonlocal loop 
operator for particle species a inserted and the RHS is obtained if we are probing the 
system in the unbroken phase governed with the groundstate denoted as and governed 
by the algebra A. We use the subscript because the value of the same diagram may be 
different if it is evaluated in a different phase with a groundstate that we will denote by 
$; in the remainder of the paper we will therefore always use brackets with a subscript. 



1.3. Topological symmetry breaking 

In this section we briefly recall topological symmetry breaking, the phenomenon that a 
phase transition to another topological phase occurs due to a Bose condensate JT7J [T_B] . 
The analogy with ordinary symmetry breaking is clear if one thinks of the particle as 
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representations of some quantum group, and assumes that a bosonic degree of freedom 
i.e. with 9 C = 1 - fundamental or composite - condenses. The breaking can then be 
analyzed, either from the quantum group (Hopf algebra) point of view, or from the dual 
or representation theory point of view |19j . 

Let us illustrate this by an example of ordinary group breaking. Suppose we have 
a gauge group 577(3) and a Higgs triplet that acquires a vacuum expectation value 
$ = (1,0,0), then the SU(2) subgroup working on the last two entries will leave $ 
invariant. Equivalently this SU(2) subgroup may be characterized by the way the 
SU (3) triplet decomposes under the SU{2) action as 3 — > 2 + 1 where the singlet 
on the right corresponds exactly to the new SU (2) invariant groundstate. In that 
sense one may select a specific residual gauge symmetry by choosing an appropriate 
Higgs representation which has a singlet under that residual group in its branching. 
For example if we want to break an SU(3) group to the SO (3) subgroup which is 
characterized by the branching rule 3 (as well as 3) — > 3 then we may choose the Higgs 
field to be in the 6-dimensional irrep of 577(3), because then 3x3 = 3 + 6—7-3x3 = 
1+3 + 5, from which follows that 6 — > 5 + 1 and again the singlet on the right corresponds 
to the SO (3) invariant vacuum state 

In the case of general quantum groups it is this branching rule approach which 
is the most natural and powerful in the context of TQFT because the fusion algebra 
corresponds to the representation ring of the quantum group. A general treatment with 
ample examples can be found in reference [19]. Let us point out some essential features 
of this procedure that one has to keep in mind. As the quantum group centralizes 
the chiral algebra in the operator algebra of a CFT, one expects that reducing the 
quantum group will correspond to enlarging the chiral algebra, and this turns out to be 
the case. In contrast to ordinary group breaking, the topological symmetry breaking 
procedure involves two steps, firstly the condensate reduces the unbroken fusion algebra 
(also called a braided modular tensor category) A to an intermediate algebra denoted by 
T. This algebra however may contain representations that braid nontrivially with the 
condensed state, i.e. with the new vacuum and if that is the case, these representation 
will be confined and will be expelled from the bulk to the boundary of the sample. 
Confinement implies that in the bulk only the unconfined sectors survive as particles 
and these are characterized by some subalgebra U C T. Let us briefly describe the two 
steps seperately. 

From A to T Assuming that a certain bosonic irrep c will condense due to some 
underlying interaction in the system, implies that c will be identified with the vacuum 
of T. For our purposes, a boson is a sector with trivial (integer) spin, though in fact 
in the context of 2 + 1 dimensions one has to also require that fusion of this field with 
itself has a channel with trivial braiding. 

The definition of the new vacuum requires to a redefinition of fields. Firstly, fields 
in A that appear in the orbit under fusion with the condensed field c are identified in 
T, so, if c x a = b then a, b — > a'. Secondly, if a field b forms a fixed point under 
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fusion with the condensate c, then the field will split at least in two parts: b — > Yli^i- 
The identifications and splittings of representations can be summarized by a rectangular 
matrix n l a that specifies the "branching" or "restriction" of fields a from in A to T with 
fields t, r,s, . . .: 



This branching matrix is a rectangular matrix (the number of particle types in the A 
and T theories is not equal in general) of positive integers. We will also consider the 
transpose of this matrix denoted as which specifies the "lift" of the fields t G T to 
fields a G A: 



One may now derive the fusion rules T from the fusion algebra ([2]). Because of 
the identifications, it is often the case that the intermediate algebra T though being a 
consistent fusion algebra, is not necessarily braided, in more technical terms, it satisfies 
the "pentagon" equation but not the "hexagon" equation. The physical interpretation 
of this fact is that the sectors in T do not yet constitute the low-energy effective theory. 
This is so because sectors t that have an ambiguous spin factor, meaning that not all 9 a 
of the lift a G t are equal, will be connected to a domain wall and hence are confined in 
the new vacuum. The confined excitations will be expelled to the edges of the system 
or have to form hadronic composites that are not confined. Yet the T algebra plays 
an important role: in [24J for example, it was shown that the T algebra governs the 
edge/interface degrees of freedom in the broken phase. 

From T tolA Some of the sectors in T will survive in the bulk, some will be confined. 
The physical mechanism behind confinement in 2 + 1 dimensional topological field 
theories is nontrivial braiding with the condensate. The vacuum state or order parameter 
should be single valued if carried adiabatically around a localized particle like excitation. 
If it is not single valued that would lead to a physical string or "domain wall" extending 
from the particle that carries a constant energy per unit length. The unconfined algebra 
U consists of the representations in T minus the confined ones, it is this algebra that 
governs the low energy effective bulk theory. The confined representations can be 
determined in the following way. First we define the "lift" of a representation in T 
as the set of representations b G A that restrict to t. Now, if all of the representations in 
the lift of t braid trivially with the lift of the vacuum, the sector t is part of U. Otherwise, 
it is confined. One may prove that the U algebra closes on itself with consistent fusion 
rules, while consistent braiding is achieved by assigning the (identical) spin factors of 
the parent sectors of the unbroken theory to the U fields. 

Let us finally mention a useful quantity, the socalled quantum embedding index q 
defined in |28j . it is a real number characterizing the topological symmetry breaking. 






a 
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This quantity is denned as 

g= E^ ; (8) 

where the index a runs over the sectors of the unbroken phase A, that correspond to lift 
of any sector u or t of the algebra U or T; the n a u is the lift of sectors u to their parents a 
and d a is the quantum dimension of the representation a. Observe that this expression 
is independent of the particular sector u, which is a non-trivial result explained in [28J. 

Choosing for u the new vacuum, we have d u = 1 and obtain that q just equals the 
total quantum dimension of the lift of the U (or T) vacuum in the unbroken A theory. 
The quantum embedding index is the analogon for the embedding index defined by 
Dynkin for the embedding of ordinary groups. As an aside we mention that the change 
in topological entanglement entropy of the disk changes also by log(D J i/D u ) = logg in 
a transition from an A to a U phase [28] . 

Let us to conclude this subsection on topological symmetry breaking illustrate the 
procedure with a very straightforward example, namely the breaking of the quantum 
group A = S , f/(2) 4 . It has 5 irreps labeled by A = 0,...,4 with spinfactors 
9 a = The A = 4 is the only boson and we assume it to condense. The 

lift of the new vacuum corresponds to the $ = + 4 of A, and hence the embedding 
index q = d + d 4 = 1 + 1 = 2. The 1 and 3 reps of A are identified, but because they 
have different spin factors, the corresponding T representation will be confined. In U we 
are therefore left with the A = 2 rep. which splits because it is a fixed point under fusion 
with the condensate as 4 x 2 = 2. We write 2 — > 2\ + 22. The values for the spin and 
the quantum dimensions and the fusion rules for these representations fully determine 
the unconfined quantum group to be U = S77(3)i. We recall that the nomenclature 
of the groups is linked to the chiral algebra, it is therefore not surprising that the 
5 , f/(2) 4 quantum group breaks to the smaller quantumgroup SU(3)i which is related 
to a larger chiral algebra. For the chiral algebras one has the conjugate embedding 
S77(2) 4 C SU(3)i which is a conformal embedding. This conformal embedding in 
turn is induced by the SO{3) C 577(3) embedding mentioned at the beginning of this 
subsection. 

1.4- Observables 

Our objective is to verify the theoretical predictions of the topological symmetry 
breaking scheme in a class of euclidean gauge theories that are expected to exhibit 
transitions between different topological phases. We will numerically evaluate the 
expectation values of various topological diagrams using Monte Carlo simulations, and 
in this section we calculate the predicted outcomes of a variety of possible measurements 
from theory. The strategy has two steps, (i) the determination of the condensate 
(including the measurement of the embedding index q) by evaluating the basic nonlocal 



open string order parameters, given by Eq. (33), (ii) measuring the socalled broken 



modular S'-matrix and from that construct the S-matrix of the U phase. We also will 
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see that the condensate fixes the branching and lift matrices and having determined 
those we can also predict the outcome of measurements of other topological diagrams 
corresponding to the lifts of U fields to A fields . 



1.4-1- Determination of the condensate and the embedding index q. We measure the 
open string operators in the model. Note that in our pictorial representation time 
flows upward, so a vertical line physically represents the creation, propagation and 
annihilation of a single particle. For the particular case of a DGT, which we study in 



this work, these lines have a realization as operators on a spacetime lattice, see Eq. (33) 
If the symmetry is unbroken we will have for any nontrivial field a that 



0. 



(9) 



$=0 



because the diagram represents the creation and subsequent annihilation of a single a- 
particle. However in the broken situation the expectation value will be nonzero for all 
fields (pi E A in the condensate which we denote by <3>. So writing, 

$ = + ^0, (10) 



we obtain that in general, 



Ocupi da ■ 



This in turn implies that it is simple to measure q as 



E 



E 



do + d ^ = Q 



(12) 



1-4-2. Determination of confinement and other topological data of the broken phase. 
Once we have determined the components of the vacuum we can determine the lifts of 
the t fields simply by studying the fusion rules of $ x a = V^i', where t! denotes the 
lifts of those t fields which contain a, i.e. for which n" = 1. Having obtained the lifts 
of the t fields the next step is to make the measurement determining whether a given t 
field is confined. This involves the measurement of the index rj, or simply: 



Ep 

aet 



E„ , 2 t if t 4 U (confined) 

^ ° ada = " dtVt = \ g*MJf t e U (not c 



confined) 



(13) 



Alternatively one can measure certain closed at loop operators that are also defined for 
fields a that split under branching and that will be defined later, for which holds that: 

if t i U (confined) 




q 2 6 u d u if t E 14 (not confined) 



(14) 



3 \"i q>j I o 

It follows that from these measurements, the fields that are confined can be determined, 
but also the quantum dimensions d u and twists 6 U of the unbroken U theory are obtained. 
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1.4-3. The broken modular S- and T -matrices. Instead of the fusion coefficients N® b an 
alternative specification of a (modular) topological field theory is by its representation 
of the modular group SL(2, Z) generated by the S and T-matrices 

S 2 = (ST) 3 = C, S* =CS = S-\ T* = T-\ C 2 = 1, (15) 

with C the charge conjugation matrix. The corresponding matrix elements can be 
expressed in the fusion coefficients and spin factors: 

s -^^^° (16) 

c 

T ab = e- 2m ^6 a 5 a , b (17) 

where T> is the total quantum dimension and the constant c is the conformal central 
charge of the corresponding conformal field theory. We recall that the central charge of 
a discrete gauge theory is zero, so in that case the T-matrix is just the diagonal matrix 
containing the spin factors. 

The great advantage of switching to the modular data is that unlike the fusion 
coefficients these generators can be directly measured using the anyon loop operators 
that arise naturally in a three dimensional euclidean formulation of the theory. We will 
evaluate the expectation value of these S-matrices numerically in our lattice formulation 
of multiparameter discrete gauge theories later on. The measured S- and T-matrix 



elements do not satisfy the relations (15) directly; however, using the measurements the 
full S- and T-matrices of the U theory, which do satisfy the modular group relations, 
can be constructed, . In the unbroken theory the measured S-matrix elements (S a b) 
correspond to the expectation values of the Hopf link with one loop colored with 
representation a and the other with representation b: 

(S ab ) = p ^ "(TJ)^ / = S ab, 

where S a b is the S-matrix of the unbroken A theory. We can however also determine the 
modular S-matrix of the residual U theory S uv directly from measurements if we take 
the splittings of certain fields a =>■ {a{\ in account appropriately. We will show how to 
do this later for the DGT's in detail and give a more general mathematical treatment 
of this elsewhere [29]. Then we will arrive at an explicit formula and algorithm to 
determine S uv : 

CLi,bj 

This expression involves not only the branching (lift) matrix n" ,; , but also the what we 
will call the broken S-matrix defined as = (5 ai (,.)$, which, because of the splitting, 
clearly involves a larger size matrix then the modular S'-matrix of the original A phase. 
From the broken S'-matrix we may directly read off S uv , the S-matrix of the effective 
low energy TQFT governed by U. An important observation is that the values of the S'- 
matrix elements in a broken phase will be different from the ones in the unbroken phase, 
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for example because of the contribution of the vacuum exchange diagram S depicted 
below, in which the condensed particle is exchanged giving a nonzero contribution in 
the broken phase while it would give a vanishing contribution in the unbroken phase: 



S, 





In the explicit calculations later on we show that this vacuum exchange diagram leads 
to a change in the S-matrix which depends on the subindices introduced above. It turns 
out that it is also possible to calculate the broken S-matrix from first principles, this 
will be discussed in a forthcoming paper [29|. 

As to be expected one finds identical rows and columns in the broken S'-matrix, 
for components that are identified, whereas the entries for confined fields will be zero. 
With this prescription the formalism outlined above is applicable in any phase of the 
theory including the unbroken one where there is no splitting and the vacuum exchange 
diagram gives a vanishing contribution. The measured T-matrix on the other hand is 
given by 



"a \KJ I 

again with (T ao ) = T a \,. After measuring or calculating the S- and T-matrices in a given 
phase, we can reconstruct the fusion coefficients with the help of the Verlinde formula 

m, 

Kb-J2 SaxS q bxS ^ - (19) 

X 

To conclude, we have in this section summarized the basic features of a TQFT and 
considered some aspects of topological phase transitions induced by a Bose condensate, 
furthermore we explained how the measurement of the L-, S-, and T-operators in 
the broken phase fully determine the quantum group of a (broken) topological phase. 
The general scheme to analyse the breaking pattern of a some multiparameter TQFT 
is to first use the open string operators to probe which fields are condensed in the 
various regions of parameter space. In a given broken phase we can subsequently 
compute/measure what we will call the broken S'-matrix S ai b., where as mentioned the 
subindex labels the splitting of the corresponding A field. From the broken S-matrix 
we can read off the S-matrix of the U theory. In the remainder of the paper we will 
explicitly execute this program for discrete gauge theories. 
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2. A euclidean lattice approach to Discrete Gauge Theories. 

2.1. Z 2 gauge theory and topological order, a prelude 

Before applying our approach to (non-Abelian) Discrete Gauge Theories (DGTs) in 
general, let us make some connections to previous work from different perspectives. 
After Wilson's seminal work on Euclidean lattice gauge theory for non-Abelian Lie 
groups to study the confinement of quarks pQ, a Hamiltonian formalism for the same 
problem was soon developed [2J. It was clear that because the gauge fields now take 
values in the group instead of the Lie algebra, one could also study models based on a 
finite group. 

These models, which we call DGTs, were mostly studied as approximations to U{\) 
or SU(N) theories in times when computers were not as powerful as today. They are 
however also interesting in their own right, since they are purely topological: there are 
no local degrees of freedom, and only the topological (generalized Aharonov-Bohm type) 
interactions survive. 

This does not automatically mean that all observables are topological quantities: 
the appearance of virtual flux-antiflux pairs gives small size-dependent corrections to 
the loop-like observables in these theories. However, since these excitations are gapped, 
these corrections are exponentially small. We will show below by explicit calculation 
that as long as one stays away from the critical points, it is justified to think of the 
observables in these theories as topological quantities. 

2.1.1. Hamiltonian formalism To connect with work other work on topologically 
ordered systems, let us first go to a Hamiltonian formalism. This is formally done 
by taking a timeslice of the spacetime lattice and taking the limit in which the temporal 
spacing goes to zero [31]. The Hamiltonian of (2+l)-dimensional Z 2 gauge theory on a 
square spatial lattice is 

H = -\xJ2( Pl - !) - E \(QpiQp2Qp3Q p a - i), (20) 

I V 

where the operators Pi and Qi act on links, the second term is a sum over the elementary 
plaquettes of the lattice where pi . . . p4 are the links of a single plaquette and A is the 
coupling constant. The operators satisfy 

{Q h P l } = 0, P? = Q* = 1, 

which means a possible representation can be given in terms of Pauli matrices Pi = cr 3 , 
Qi = <j\ acting on spin-^ bosons living on the links. Note that the algebra above is 
the same as the Z 2 version of ([I]), and indeed a closed string of Pj operators generates 
a Wilson loop, whereas a closed string of Qi operators creates a closed Dirac-string. 
Gauge transformations act on the star of four links il . . . i4 adjacent to a site i 

Gi = PnPi2PaPni 
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and build the gauge invariant Hilbert space, one has to implement a Gauss law for 
physical states 

(1-Gi)|^) = for all sites %. (21) 

Now we can make the connection with work by Kitaev [32] and Wen [33] . Their models 
(Toric code, Z 2 string nets) correspond to Hamiltonian Z 2 DGT where the coupling 



A = and the gauge constraint (21) is not strictly enforced. Setting A = makes the 
theory purely topological, the ground state is an equal weight superposition of all states 

n>i°>> ( 22 ) 

lac 

where C is a closed loop of links and |0) is the state with the property Pi\0) = |0) for 
all links I. Viewing the link variables as spin-| bosons, this vacuum state corresponds 
to all the spins being in the up state. Since the expectation value of any loop operator 



(22) in the ground state is equal to one, and these loops are Wilson loops in the gauge 



theory language, the theory is topological. 



By not enforcing the gauge constraint (21) strictly but adding it as a term to the 



Hamiltonian, these models allow for massive open strings. Such open strings are not 
gauge-invariant at their endpoints and therefore correspond to external charges. 

2.1.2. Euclidean formalism The Z 2 gauge theory in the Euclidean approach, where we 
discretize both space and time, is described by the action 

S = -PJ2 U p^ U p^ (23) 

v 

where the sum is again over all plaquettes (now both spatial and temporal) and the U 
variables are numbers ±1. The partition sum 

z = -£'-* 

{U} 

and the expectation value of gauge invariant operators O 
{u} 

are the quantities of interest here. The gauge invariance, which in the Hamiltonian 
formulation was enforced by projecting out states from the Hilbert space, is now manifest 
in the action and the operators. The partition sum is over all gauge field configurations, 
but since all sums are finite, gauge fixing is not required^] 

If the coupling f3 is large, the dominant contribution from the partition sum will be 
from field configurations where all plaquettes UUUU = +1. In the limit /? — > oo this 
is strictly true, and one is left with a topological quantum field theory, as was the case 



for the Hamiltonian (20) with A = 0. For (3 small there is a confining phase, the phase 



transition is at f3 = 0.7613 [34J. 

| This even holds for continuous groups, since we integrate over the group instead of the algebra. 
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In most of this work, we study the topological properties of a DGT, for a general 
group H. To show that for finite coupling constant (3 this is good approximation, let us 
perturbatively calculate the expectation value of a Wilson loop in this Z 2 theory. The 
Wilson loop W{C) is the product of U variables around a closed loop C 



(w(c)) = ±J2 uu --- Ue ~ 



z 

{u} 

For large /3, the action is minimized by configurations for wich all plaquettes are +1. The 
first order perturbation comes from those configurations in which one link is —1. In three 
dimensions, this excites 4 plaquettes, so the Boltzmann weight for such configurations 
is e~ 4/3 smaller than for those with no excited plaquettes. 

If the lattice has size N x N x N, there are 3N 3 links. For a contour C of length 

This shows the corrections to the purely topological result W{C) = 1 are, for (3 several 
times larger than the critical point, negligible for simulations of reasonable lattice sizes: 
for a Wilson loop size 10 x 10, (3 = 3.0 yields corrections only in the third digit. 

Another gauge-invariant quantity is the 't Hooft loop, which lives on a loop C of 
the dual lattice. Such a loop pierces a number of plaquettes p, and the 't Hooft operator 

H{C) = J^J e~ 2/3 u pi u p2 u p3 u p* 
pec" 

flips the sign of the coupling (3 —> —(3 for these plaquettes. This forces a Z 2 magnetic 
flux through these plaquettes. We will define operators generalizing the 't Hooft and 
Wilson loops for general non-Abelian DGTs shortly. 



2.1.3. Phase structure The action (23) can realize three phases when one also allows 
for negative coupling. For large positive /3, the phase mentioned before is realized, where 
almost all plaquettes are +1. For large negative (3, almost all plaquettes are —1. For 
small |/3 1, a confining phase where the magnetic Z 2 flux has condensed is realized. 

To study the phase diagram of a (non-Abelian) DGT in full, we find it convenient 
to formulate the action in the class basis, instead of the irrep basis. This means we 
do not take the character of the group element of the plaquette product UUUU in the 
action, but we define delta functions on each class. We will explain in detail how this 



works in section 2.3 For Z 2 the phase diagram is one-dimensional, but the introduction 



of a second coupling constant will get rid of the need for negative couplings: 

S = - (P+iS+i(U P iU p2 U p3 U p4 ) + (3^5^(U pl U p2 U p3 U p4 )) , 
p 

where Sa{U) for a group element U and a conjugacy class A gives +1 if U G A and 
zero otherwise. For non-Abelian groups this formulation makes the phase diagram much 
more intuitive, for Z 2 it is rather artificial. In Figure [T] the phase diagram of the pure 
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Figure 1. Phase diagram for a pure Z2 gauge theory. 

Z 2 gauge theory is shown as a function of the conjugacy class couplings (3 + i and 
Later in this work we present similar phase diagrams for the D2 gauge theory. 

It is well-known that the inclusion of matter coupled to the gauge fields complicates 
the phase diagram strongly. The question of whether there exist good order parameters 
to distinguish the phases in coupled gauge-matter systems is interesting in its own right 
and highly non-trivial [33], but it is not something we will go in to here. 

2.2. DGT and the quantum double of a finite group 

The particles in a Discrete Gauge Theory, their fusion and braiding properties, spins 
and so forth are all obtained by working out the representation theory of the underlying 
quantum group, which is the quantum double of the finite discrete subgroup. We will 
not give a detailed account on the emergence of quantum group symmetry in DGT, this 
can be found in the literature |14| . but do present a short summary of the basics to fix 
the notation and introduce some key concepts required later on. 

Consider the following operators acting on states in the Hilbert space of a DGT. 
First there is the flux projection operator, denoted by Ph, which acts on a state 



Secondly, we have the operator g, for each group element g G H, which realizes a global 
gauge transformation by the element g: 



where it should be noted that we have not yet modded out by the gauge group to obtain 
the physical Hilbert space. These operators do not commute, and realize the algebra 




if the state \ip) contains flux h 
otherwise 



g\ip) = \ g i>) 



PhPh' 
gP h 



— &h,h'Ph 

= Pghg-^9 ■ 
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The set of combined flux projections and gauge transformations {Ph.g}h, g &H generates 
the quantum double D(H), which is a particular type of algebra called a Hopf algebra. 

The representation theory of the quantum double D(H) of a finite group H was 
first worked out in [15] but here we follow the discussion presented in |36j and follow 
the conventions of those lecture notes. 

Let A be a conjugacy class in H. We will label the elements within A as 

for a conjugacy class A of order k. In general, the centralizers for the different group 
elements within a conjugacy class are different, but they are isomorphic to one another. 
Let A N C H be the centralizer for the first group element in the conjugacy class A, 
denoted by A h\. 

The set A X relates the different group elements within a conjugacy class to the first: 

A X = { A x ftl ,4c h2 , . . . ,^ hfe | % = % % A x~ h l) . (24) 

This still leaves a lot of freedom, but we fix our convention such that A Xh 1 = e, the 
group identity element. The centralizer A N, being a group, will have different irreps, 
which we label by a. The vector space for a representation a is spanned by a basis 
a Vj. The internal Hilbert space corresponding to an irrep of the quantum double that 
combines magnetic and electric degrees of freedom, V( A ' a \ is then spanned by the set 
of vectors 

{|V«,-}} , 

where i runs over the elements of the conjugacy class, i = 1,2, ... , dim A and j runs 
over the basis vectors of the carrier space of a, j = 1, 2, . . . ,dim a. These irreducible 



representations correspond precisely to the particle types a,b, . . . in section 1.2 They 
obey a set of fusion rules as in ^ and it is possible to calculate the modular S'-matrix, 
F-symbols and so on. 

To see that this basis is a natural one to act on with our flux measurements and 
gauge transformations, consider the action of a pure flux projection P^e 

P h e\%, n Vj) = 8^%,^) , 

and a pure gauge transformation J2h ^g 

^PhglXfv,) = \gXg~\Y, D ^9)mj a v m ). 

h m 

The matrix action ]T A ' a ) of an irreducible representation (A, a) of some combined 
projection and gauge transformation P^g: 

rt A ' a \P h g)\X, a Vj ) = 6 h>g A hig -i\gXg-\J2 D °<&™i av ^ ' ( 25 ) 

m 

where the element g is the part of the gauge transformation g that commutes with the 
flux A hi, defined as 

~9=\l 9 -^ A x h ^ (26) 



The modular S -matrix as order parameter for topological phase transitions 



18 



's 




Figure 2. In our convention a gauge transformation g at location i\ transforms — > glJi^, 



and D a (-)ij is the matrix representation of the centralizer. 

To conclude we give a simple expression for the modular S'-matrix that can be 
obtained by calculating the trace of the monodromy matrix 

S( A ,a){B,p) = TJf\ Traix^hxg)* Tr ^(xl 1 gx h )* , (27) 

g&A,h£B,[g,h]=e 

where [g, h] is the group theoretical commutator: [g, h] = ghg~ x h~ l . 
2.3. Lattice actions and observables 

We discretize three-dimensional spacetime into a set of sites i, j, ■ ■ ■ using a rectangular 
lattice. The gauge field Uij, which takes values in the gauge group H, lives on the links 
ij, jk, . . . connecting sets of neighboring sites. The links are oriented in the sense that 

'V cp- 

We note that the gauge field Uij takes care of the parallel transport of matter fields 
that are charged under the gauge group from site i to site j. An ordered product of 
links along a closed loop is gauge invariant up to conjugation by a group element and 
measures the holonomy of the gauge connection. Gauge transformations are labeled by 
a group element ^ G H and are performed at the sites of the lattice. The gauge field 
transforms as 

Uij ^ giU i:i gj 1 , (28) 

where the orientation of the links (incoming or outgoing) has to be taken into account 
as shown in Figure [2] 

The standard form for the lattice gauge field action makes use of the ordered product 
of links around a plaquette ijkl: 

Up = Ujki = U^ Ujk Uki Uu , 

which transforms under conjugation by the gauge group, 

U p ^ gi.U p gi l . 
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The gauge action per plaquette which corresponds to the Yang-Mills form F£ v in the 
continuum limit for H = SU (N), is given by 

S P = -J2^Xa(U p ) , (29) 

a 

where Xa is the group character in irrep a and (3 a is inversely proportional to the square 
of the coupling constant for irrep a. This is known, for H = SU (2) and the sum over 
representations limited to the fundamental one, as the Wilson action pQ. The action 



(29) is the euclidean analog of the Kitaev Hamiltonian introduced in [32] (which in 



itself is a variation of the Kogut-Susskind Hamiltonian [2]) where the electric field is 



represented by the timelike plaquettes and the gauge constraint (21) is not necessary 
since everything is manifestly gauge invariant. 

For SU(N) gauge theories one usually only includes the fundamental representation 
and is thus left with only one coupling constant. This is not necessary however: 
gauge invariance of the action is ensured by the fact the characters are conjugacy class 
functions, in and therefore we will consider actions where the number of independent 
couplings equals the number of conjugacy classes i.e. the number of irreps (for a finite 
group these numbers are finite and equal). 

For our purposes, namely the study of magnetic condensates in DGTs, equation 



(29) is not the most convenient to work with. We perform a change of basis in the space 
of coupling constants to write it as a sum over delta functions on conjugacy classes: 
SA{h) — 1 if h G A, and otherwise. In this basis the action becomes 

A 

This formulation allows us in particular to directly control the mass of the different 
fluxes in the theory, which will ease the search for different vacua in the phase diagram. 
Increasing the coupling constant for a certain conjugacy class (magnetic flux) A will 
increase the contribution of configurations carrying many A fluxes to the path integral. 
Likewise, setting all (3 a to zero except (3 e , the coupling constant for the trivial conjugacy 
class, will result in an "empty" vacuum and therefore an unbroken phase. 

To perform the transformation to the conjugacy class basis, we need to make use 
of the following orthogonality relations valid for all finite groups H 

dg Xa(g)x*p(g) = <W, (30) 

H 

Y,^{g)x* a {h) = ^Sg,heA (31) 
= otherwise, 

where \H\ is the order of the group H, \A\ is the order of the conjugacy class A, 1Z is 
the set of irreps and group integration is defined as 



l H dgf{g) = ^Y,f^)- 
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Equations (30) and (31 ) show that the irreducible representations of a group H form an 
orthonormal set for functions on conjugacy classes of H. We thus expect the conjugacy 
class delta function to be expressible in terms of characters 

for some set of constants {c a }. We multiply both sides of this expression by a character 
of the same group element in another irrep (5 and perform the integrations by use of the 
orthogonality relations d30| and (J3ij) 



n 



x*M) 



W\ 



C a S a /3 — C/3, 



dg x}(g)xa(g), 



where the slightly abusive notation Xa{A) means the character of any group element of 
A in the representation a. This shows that 



A\ 



(32) 



which in turn implies that the the difference between (29) and (2.3) is just a change of 
basis: 

Aec a^n 
where C is the set of conjugacy classes and is given by 

a 

To probe the physics of the system for a fixed set of values of the coupling constants 
in the action, we will use a set of order parameters and phase indicators. These order 
parameters are in one-to-one with the set of fundamental anyonic excitations of the 
theory. 



Order parameters and phase indicators. We distinguish two different sets of order 
parameters that are closely related to one another. The first is the set of closed loop 
operators, that physically correspond to the creation, propagation and annihilation of 
an anyon-anti-anyon pair in spacetime. The second is the set of open string operators 
that create, propagate and annihilate a single anyon. In the background of a trivial 
vacuum, only the loops can have nonzero expectation values, since the creation of a 
single particle would violate the conservation of the quantum numbers of the vacuum in 
such a background. This means that the open strings tell us something about possible 
Bose condensates, whereas the closed loops tell us about the behaviour of external 
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Figure 3. The first two plaquettes appearing in expression (33). The ordered product U p of links 
around a plaquette p needs to be taken with an orientation that has to be constant throughout 
the loop. 



particles put into this background. We define the open string operators only for the 
purely magnetic sectors, since in this work we only study magnetic condensate^} 

First we will define the loop operators. This set of nonlocal order parameters was 
introduced in a previous publication [16J. For a full discussion, we refer to that work. 
Here we recall the essentials and fix the notation. The closed loops are a generalization of 
the Wilson and 't Hooft loops. They create a particle-antiparticle pair from the vacuum 
and annihilate them at a later time. These loops allow us to calculate Aharonov-Bohm 
type phases and determine which anyonic excitations will be confined. In SU(N) gauge 
theories, the Wilson loop for a free excitation, e.g. in the Higgs phase of SU (2) theory, in 
general falls off as e~ cP , with P the perimeter of the loop, whereas a confined excitation, 
such as the 3 charge of an external quark source in pure SU (3) gauge theory describing 
QCD, falls off as e~ c ' A , with A the area of the loop. 

Because the excitations in a DGT are gapped, numerically we find that the 
expectation values of loop operators are constant as a function of size. The argument 



for this behaviour for the Z2 theory is in Section 2.1.2 Although only strictly true in the 
limit of infinite coupling constant, the gap supresses the dependence on size so strongly, 
we will assume that the theory is a purely topological one in the region of coupling 
constant space in which we are interested. 

Let us draw a closed loop on the dual lattice, this loop pierces a set of plaquettes 
C through which we will force magnetic flux. Now draw another loop , this time on 
the real lattice, such that (i) each point of this loop lies on the corner of a plaquette in 
C and (ii) the two loops do not linhjJQ The combination of the two loops establishes a 

§ Electric condensates break the gauge group H to some subgroup K by the conventional Higgs 
breaking, this implies in the present context that D(H) will be broken to D(K), which in turn means 
that the fluxes in the coset G/H are confined 18J. 

|| One can also create loops that have a linking number, the expectation values of such loops allow one 
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Figure 4. A set of plaquettes forming a closed loop on the lattice. The fat links constitute the 
/i-forest. 

framing: we have selected a location for the electric charge of a flux-charge composite. 
This framing also provides us with a point and orientation on each plaquette from which 
to take the plaquette product U p : for non-Abelian groups the product of the four links 
depends on which corner you start. 

To insert a flux h in a plaquette p G C, we have to "twist" the Boltzmann factor of 
this particular plaquette by locally changing the action from S(U P ) to S(h~ 1 U p ): if the 
minimum of the action was previously obtained for U p = e, it is now shifted to U p = h. 
We want to perform this twisting procedure for all plaquettes in C. 

The notion of a group element as a magnetic flux is not gauge- invariant: under a 
gauge transformation by g, a flux h transforms as ghg~ x . Therefore it is necessary to 
sum over the group elements h in a conjugacy class A in some way. 

One can go about this in two different, and inequivalent, ways. 

• The authors in Ref. [37J only studied pure magnetic flux loops (without electric 
charge) and performed a sum over the conjugacy class for each plaquette in C 
individually. This leaves a gauge-invariant expression, but the loop loses its framing, 
since a conjugation by the element {7 i maps 

h- 1 U 01 U 12 U 23 U 30 -> U 01 h- 1 U 12 U 23 U 30 , 

for a plaquette spanned by group elements Uqi . . . U 3 q. 

• When dealing with nontrivial braiding properties of loop operators it is necessary 
to choose a basepoint io in space with respect to which all operators are defined, it 
provides a calibration that serves as a "flux bureau of standards" , borrowing a term 
from [38J. This point can be anywhere in spacetime and does not need to be on 
the loop. We then define a function k ip (h, {Uy}) of h and the gauge field variables 
{Uij} for the twist element that has to be inserted into the plaquette product for 

to calculate the topological spin for a given excitation. 
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plaquette p, where i p is the corner of the plaquette chosen in the framing 

K = K(h,{u ij }) = ur o ih- 1 u ioip . 

With this notation and the above considerations, the anyonic operator is 
given by [J 

A<^»(C) = J2U D " Ul^vA, "I-* **) (33) 

Here iterates over the plaquettes in C and D a is the representation function of 
the centralizer irrep a of A N. The link Uj-ij neighbours the plaquette pj, and the 
combination in brackets always takes values in the centralizer subgroup of the conjugacy 
class A. The exponential of the difference of two actions changes the minimal action 
configuration to one containing flux h for the plaquette under consideration. 



The operator in expression (33) is a generalization of the Wilson and 't Hooft 
loops, and by constructing it we have established the desired one to one correspondence 
between irreducible representations of the quantumgroup and loop operators for the 
pure discrete gauge theory. If we fill in for A the trivial conjugacy class, the exponent 
vanishes and the x group elements are equal to the group unit, so after we multiply out 
the D Q -matrices we are left with 

where the product of Us is an ordered product along the loop on the lattice. 

On the other hand, if we replace a by the trivial representation, we are left with 

A^\C) = II e s{u ^- s(hPWp i\ 

heApj&C 

which is comparable to the order parameter proposed in [37], but the gauge invariance 



with respect to the transformations (28) is ensured in a different way. We sum over 
the conjugacy class only once and insert the flux in a gauge invariant way by parallel 
transporting it along the loop from a fixed basepoint. The operator in [37] sums over 
the conjugacy class for each individual plaquette. This way also gauge invariance is 
achieved, but the loop loses its framing, and therefore is not suitable to describe true 
anyonic charges. 



The open magnetic string operators are a variant of expression (33) where the set 
of plaquettes C corresponds to an open string on the dual lattice. Looking at the h- 
forest configurations, it can immediately be seen that such a string, corresponding to the 
creation and subsequent annihilation of a single particle, has zero expectation value in 
the trivial vacuum. For these strings to acquire a non zero expectation value a vacuum 
exchange contribution is required, which we will focus on now. 

% This definition is different from our original definition 16 by a factor of pr. This definition gives 
the correct S'-matrix elements directly. 
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Figure 5. The vacuum exchange contribution. Double /i-forest, configuration contributing 
to the S'-matrix measurement of two non commuting fluxes. The fat links are the g-forest and 
/i-forest and the shaded plaquettes are a string of [g, h] flux connecting the the two loops. 



The vacuum exchange contribution. We use the set of operators {A^' a ^} to measure 
the elements of the S'-matrix by picking two loops C\ and C<i that link each other once 

(S { A,a)(B,p)) = (A^C^A^)^)) . (34) 

In the trivial vacuum the 5 , (A,a)(s,/3)-matrix elements of fluxes g G A and h G B for which 
ghg' 1 h' 1 = [g, h] ^ e evaluate to zero (this is what we measure using the operators 
(33) and calculate algebraically (27) ). If we however measure the S'-matrix elements of 
such noncommuting fluxes in a broken vacuum nonzero matrix elements can appear. 

This is most easily explained by considering an example. The main contribution to 
a single loop of pure magnetic flux is of the form pictured in Figure |4j This configuration 
is called the h-iorest state in earlier literature [37]. Modulo gauge transformations this 
is the dominant configuration in the trivial vacuum that contributes to a loop of flux 
labeled by conjugacy class A, where g G A. Expression (33) contains a sum over these 
group elements within a conjugacy class, but let us for now focus on one of the group 
elements. Each link in this configuration has value e, except for the fat links in Figure 
|4j they have value h. That this configuration leads to a loop or tube of flux is easily 
seen: within the forest each plaquette has a value eh eh' 1 = e, whereas at the edges 
the value is ehee = h (depending on the orientation of the plaquette product). This 
is also the easiest way to see the origin of the Aharonov-Bohm effect on the lattice: an 
electric charge loop having linking number 1 with the flux loop will have exactly one 
link with value h in it, therefore its value will be Xa(h)- 

Consider now the dominant configuration that contributes to the S'-matrix element 
S(a,i){b,i)- We again pick two group elements g G A, h G B and draw a similar diagram. 
This is shown in Figure [5j By similar logic this causes the plaquettes at the boundary 
of either forest to have value g respectively h. Inside the forests most plaquettes still 
have value e, however there are some plaquettes that are different. There is a tube of 
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plaquettes that have value [g,h], where the two forests intersect. In general this group 
theoretical commutator is not equal to identity element for nonabelian groups. This 
is the physical reason behind the appearance of zeroes in the S-matrix for nonabelian 
theories. This tube of plaquettes represents a flux [g, h] going from the one loop to the 
other. In the trivial vacuum this flux will be gapped, so the contribution of this diagram 
to the path integral expectation value will be negligible. 

However, a different situation appears when we are in a vacuum where the flux [g, h] 
has Bose condensed. We cannot give a single configurations that contributes dominantly 
to the path integral (there are many), but we can say that configurations like the one 
in Figure [5] are now contributing since the mass for the flux [g, h] has disappeared. 

Thus we expect that in the measurements there will be cases where zeroes in the 
original S'-matrix will obtain a nonzero value in the broken phase. 



An auxiliary A N gauge symmetry. The operators (33) are invariant with respect to the 
local H gauge transformations (28). However, in our formulation of the operators we 
have tacitly introduced another, auxiliary A N gauge symmetry that is less obvious. A 
crucial property that allows one to determine the topological symmetry breaking pattern 
in detail is that the loop operators do transform nontrivially under this symmetry. In 
a non trivial ground state, these symmetries may be broken and will therefore lead to 
the lifting of certain degeracies related with the splitting of fields in the topological 
symmetry breaking process. So this hidden symmetry turns out to be a blessing in 
disguise. 

Let us first note that there is no preferred choice for the coordinate system (24) we 
define for the conjugacy classes. Once a certain choice {x^} has been made such that 
hi = x^hix^, a set {x' h .} with 

x 'hi = %h t n hl , [n hi , hi] = e, (n hi e A N) (35) 

will do just as well. In the trivial vacuum, the S'-matrix is invariant with respect to this 
transformation. This is most easily seen by looking at the algebraic expression (27), but 
it is also confirmed by our measurements of (34). 

This invariance can be understood on the operator level by multiplying out the 
representation matrices of the centralizer in equation (33). Generally this will lead to 
terms of the form 

Tr a <7 = Tr a (x^gx hi ), 

where g is the product of links on the loop and h\~ = ghtg -1 , implying that indeed 
g G A N . When the loop is linked with another loop, the element g will in general be in 
the conjugacy class of the flux of this other loop. Under the transformation (35) of the 
conjugacy class coordinate system, the above expression will transform as 

Tl a(n^x^gx ht n hi ) = Tr a (n h .n^x^gx hi ), 

due to the cyclicity of the trace. This elucidates the invariance of the S-matrix in 
the trivial vacuum under the translation of the x^:. non-commuting fluxes never have 
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a non-zero matrix element, and if [g, h] = e, we have that hi = hf. and therefore 
n h l n ~hl = e - 111 a non-trivial ground state where non-commuting fluxes may have non- 
zero S'-matrix elements due to a vacuum interchange contribution, the transformation 



(35) may manifest itself in different measured matrix elements. This means that in such 



cases the entry (A, a) may split into multiple entries {(A i: As we are interested in 
these multiple entries, we will in our calculations always include the nontrivial behavior 
of our observables under this auxiliary A N action. This turns out to be one of two 
mechanism responsible for the splitting of irreps of A into multiple irreps of U, the 
other of which we turn to now. 



An auxiliary H/ A N symmetry. There is another symmetry, but now on the level of 
the fusion algebra that turns out to be useful. Suppose in the theory there exists a rule 
of the form 



(A,a) x (e,/3) = (A, a), 



(36) 



where (e, /3) is some one-dimensional purely electric representation. This turns out to 
be the case whenever the representation IT- 6 '^-) evaluates to unity for all elements in 
A N, the normalizer of conjugacy class A. We can prove this using the explicit expression 
for the fusion coefficients in terms of the quantum double characters: 



N (A,a)(BJ3) 
iV (C,7) 



-L Tr [ n(Aa) ® n(B,/3) ( A ( p h g))] Tr [n(^) (p h g 

g,h 



(37) 



Picking (B,fi) = (e,p) and ((7,7) = {A, a), 



N, 



(A,a)(e,/3) 
(A,a) 



u (A,a) ^ n (e,/3) 



2^ P hl g®P h2 g 

\h\h2=h 



g,h 

■■■TT[U^\P h g)]* 
]- VTr [rt A > a \P h g)®IlW\P e g)] Tr [U^ a \P h g 



\H\ ^ 

g,h 

1 

\H\ 



J2 Tr [U^ (P h g)] Tr [U^ (P h g) 

g& A N,h&A 

where in the latter line we have made use of the orthogonality of the characters. We 
assumed Ii^\P e g) = 1 for all g e A N. The sum over h is restricted since if h ^ A the 
matrix element will be zero and the sum over g is restricted since if g qL A N the matrix 
element will be off-diagonal and thus not contribute to the trace. 



So we see that the fusion rule (36) leads a degeneracy in the calculation of S'-matrix 
elements since by definition 




(A,a) (CVy) 
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However, on the operator level this equality does not hold. Indeed, when we probe the 
LHS of this equation in a non-trivial vacuum the result will in general differ from the 
RHS. In particular, it turns out that the different U representations that lift to the 
same A representations (A, a) differ precisely by such a fusion. So this degeneracy may 
be lifted in the broken phase and give rise to a additional splittings of certain entries 
(A, a). Consequently in our numerical calculations we have to explicitly keep track of 
the presence of such electric representations (e,/3), that satisfy (36) and see whether 
they give rise to additional splittings. 

To conclude this section, we remark that we have very explicitly indicated how one 
gets from the modular S- matrix S a b to the extended or broken S"-matrix 5 0i &., from 
which the topological data of the broken U phase can be immediately read off. 



3. The D(D 2 ) gauge theory 

We turn to the particular example we have chosen to work out in detail: a discrete gauge 
theory with gauge group D 2 , also called the quaternion group. The representation theory 
was worked out in [36J, here we summarize the results that are required to describe the 
breaking by a Bose condensate. 



3.1. Algebraic analysis 

The group D 2 contains eight elements that can be represented by the set of 2 x 2-matrices 
{t,—t,±iai,±ia 2 ,±io- 3 }, (38) 



where the ov 



1, 2, 3 are the Pauli spin matrices. We denote the conjugacy classes 



as e = {1}, e = {-1}, X x = {ia u -ici}, X 2 = {i(r 2 ,-ia 2 }, X 3 = {ia 3 , -ia 3 } 
and the irreducible group representations as 1, the trivial irrep, Ji, J 2 , J 3 three one- 



dimensional irreps and x the two-dimensional irrep given by (38). The character table 



is given on the left hand side of Table [T] The centralizer groups for the conjugacy 



D 2 


e 


e 


x l 


x 2 


x 3 


I 


I 


I 


I 


I 


1 


■h 


1 


I 


I 


-I 


-1 


h 


1 


I 


-I 


I 


-1 


h 


1 


I 


-I 


-I 


1 


X 


2 


-2 












z 4 


1 


ioi 


-1 


-ioi 


r° 


1 


I 


I 


I 


r 1 


1 


i 


-I 


—i 


r 2 


1 


-1 


I 


-1 


r 3 


1 


—i 


-I 


i 



Table 1. Character table of the group D2 and of Z4 as a centralizer of the conjugacy class Xj. 

classes e and e are both D 2 since the elements in these conjugacy classes constitute the 
center of the group. The conjugacy classes Xj, i — 1, 2, 3 have non-trivial Z 4 centralizer 
subgroups, of which the character table is given on the right hand side of Table [Tj The 
irreducible representations of the quantum double are labeled by a combination (A, a) 
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of a conjugacy class A and a centralizer irrep a. The full set of fusion rules for the 
D{D 2 ) theory is given in Appendix A All in all, there are 22 sectors: the trivial flux 
paired with the five irreps of _D 2 , the e flux paired with the five irreps of _D 2 and the 
three Xj fluxes paired with the four Z4 irreps. The sectors that involve an Xj flux or 
a x irrep have quantum dimension 2, the others have unit quantum dimension. One 
obtains that the total quantum dimension for the theory Da = 8. 



Breaking: (e, 1) condensate. In this case the lift of the new vacuum is = (e, 1) + 
(e, 1), which implies that q = d( e ,i) + rf(§,i) = 2. To determine the effective low energy 
theory we fuse <fi with all particle sectors of the theory and look for the irreducible 
combinations that appear. As before the notation (A, a) stands for a particle with 
magnetic flux A and electric charge a. 



4> X ( e ; 1) 


= (e,l) + (e,l) 


(j) x (e, Ji) 


= (e, Ji) + (e, Ji) 


<P x (e,x) 


= (e,x) + (e,x) (*) 


4> x (e, 1) 


= (e,l) + (e,l) 


<P x (e, Ji) 


= (e, Jj) + (e, J) 


<P x (e, X) 


= (e,x) + (e,x) (*) 


x (x 4 ,r ) 


= (x J ,r°) + (x i ,r ) 


x (A^r 1 ) 


= (x i ,r 1 ) + (A i ,r 3 ) (*) 


x (X t ,T 2 ) 


= (x i; r 2 ) + (x 4 ,r 2 ) 


x (X U T 3 ) 


= (x i ,r 1 ) + (A i) r 3 ) (*) 



The lines marked with (*) have components on the right hand side that carry different 
spin factors, implying that they are confinement in the broken phase. Studying the 
fusion rules of the surviving combinations of irreps leads to the conclusion that the 
effective U theory is Z?(Z 2 <S> Z 2 ). We denote the four different irreps and conjugacy 

classes of the group Z 2 <8> Z 2 by the labels ++, H — , — h, , the first (second) symbol 

standing for the first (second) Z 2 . This means T>% = 32 and T>^ = 16. The branchings 
of A irreps into the unconfined U theory are 



(e, 1) + (e, 1) 




(++> 


+ + ) > 


^(++,++) 


= 1 


e, Ji) + (e, Ji) 


— > 


(++, 


+ -) > 




= 1 


e, J 2 ) + (e, J 2 ) 


-> 
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-+) > 
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e, J 3 ) + (e, J 3 ) 


— > 


(++> 




) 


= 1 




— > 


(-+> 


++) > 




= 1 


(A^r ), 


—> 


(-+> 


+ -) , 




= 1 


(A^T 2 )! 


-> 


(-+: 


- + ) > 




= 1 


(A 1; r 2 ) 2 


-> 


(- + > 




d i-+,—) 


= 1 


(a 2 , r°)x 


-> 


( + -> 


++) > 




= 1 
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(x 2 


r u ) 2 


->■ 


(h — 


-+) , 


7 

«(+- 


-+) 
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(x 2 


r 2 )i 


->■ 


(h — 
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.+-) 
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(x 2 


r 2 ) 2 


->■ 


(h — 


— ) , 


7 

«(+- 


,— ) 


= 1 


(X 3 


r )! 


->■ 




++), 




,++) 
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(x 3 


r°) 2 


->■ 




— ), 




,— ) 


= 1 


(x 3 


r 2 )i 


->■ 




+-), 


d( 


,+-) 


= 1 


(X 3 


r 2 ) 2 


->■ 




-+), 


d( 


-+) 


= 1 



which all have quantum d 



mension cL 



(e, X) + (e, x) -> h , dt! 



d 



(x 1 ,r 1 ) + (x 1 ,r 3 ) t 2 , 
(x 2 ,r 1 ) + (x 2 ,r 3 ) * 3 , 
(x 3 ,r 1 ) + (x 3 ,r 3 ) t 4 , d u 



d t3 = 2 



1 , while the confined fields are 

2 
2 



and have d t — 2. 



Breaking: (Xi,T°) condensate. There is an obvious symmetry in the fusion rules 
between the three (X iy T°) particle sectors. We choose to study the case where the 
(Xi, r°) condenses. This gives for the new vacuum = (e, 1) + (e, 1) + (X 1? r°), from 
which follows that q = 4 in this case. We now read off the lifts of the T fields on the 
right: 



x (e, 1) 


= (e,l) + (e,l) + (X 1 ,r°) 


x (e, Ji) 


= (e, Jj) + (e, Jj) + ^(Xx, r°) + r hi (X 1 , T 2 ) 


x (e,x) 


= (e,x) + (e,x) + (x 1 ,r 1 ) + (x 1 ,r 3 ) 


x (e, 1) 


= (e,l) + (e,l) + (X 1 ,r°) 


x (e, Jj) 


= (e, Ji) + (e, Jj) + 8 U {X U r°) + 77h(^i, T 2 ) 


x (e, x) 


= (e,x) + (e,x) + (x 1 ,r 1 ) + (x 1 ,r 3 ) 


0x (x 1; r°) 


= (Xx, r°) + (X u r°) + (e, 1) + (e, 1) + (e, J,) + (e, J x ) 


0x (x^r 1 ) 


= (X 1 ,r 1 ) + (X 1 ,r 3 ) + (e,x) + (e,x) 


0x (X!,r 2 ) 


= (X u T 2 ) + (X u T 2 ) + (e, J 2 ) + (e, J 2 ) + (e, J 3 ) + (e, J 3 ) 


0x (X!,r 3 ) 


= (X 1 ,r 1 ) + (X 1 ,r 3 ) + (e,x) + (e,x) 


0x (x,,r°) 


= (x,r°) + (x l ,r ) + (x fc ,r°) + (x fc ,r 2 ) (i^fc^i) 


0x (x^r 1 ) 


= (x, r 1 ) + (x, r 3 ) + (x k , r 1 ) + (x fc , r 3 ) 


0x (x,,r 2 ) 


= (x^ r 2 ) + (x^ r 2 ) + (x fc , r°) + (x k , r 2 ) 


0x (x,,r 3 ) 


= (x, r 1 ) + (x^ r 3 ) + (x k , r 1 ) + (x fc , r 3 ) 



We have used the symbol 5ij which is 1 when i and j are equal and is zero otherwise, 
and rjij which is 1 when % and j are not equal and is zero when % and j are. The U theory 
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is D{Z 2 ) 
fields are: 



Z 2 ® Zo. This means T> 



T 



16 and 



(e,l) + (e,l) + (X 1; r c 



(e, Ji) + (e, Ji, 



(*i,r°) 2 
(x 2) r ) 1 + (x S) r°) 1 
(x^r^ + p^r 2 )! 

and of the confined fields: 
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4. The lifts of the unconfined 
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— >• 


, d tl 
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= 1 




) + (x 3 ,r 3 ) 


->• 


t& , d tfi 


= 2 



3.2. Measurements by lattice Monte Carlo simulations 

The five couplings {/3a} for conjugacy class A that appear in the action of the D(D 2 ) 
theory 

S p = £ - {/WOp) + /3e5e(C/ P ) + ^x 1 S Xl (U p ) + Px 2 6x 2 (U p ) + (3x 3 8x 3 (U p )} , (39) 
p 

are inversely proportional to the masses of the fluxes A. For example if we put all 
couplings to zero except for (3 e , which we make large (at least as large as 2.0 as we 
will see shortly), the trivial vacuum is realized: this is the configuration where for 
all plaquettes U p = e. Deviations from this configuration occur because of quantum 
fluctuations, but since all excitations are gapped they will be exponentially suppressed. 
The gap in this vacuum is easily calculated to be of the order of 4/3 e , since the smallest 
excitation above the configuration in which all plaquettes are e is one in which one link 



has a value ft / e. This excites four plaquettes and changes the action (39) by a value 
of 4&. 



Monte Carlo considerations. For the other, nontrivial phases in this theory, the 
dominant configurations contributing to the path integral are not so readily identified. 
To gain insight into what configurations contribute we use a Monte Carlo simulation, 
in particular a modified heat bath algorithm. Bluntly applying this algorithm to our 
problem leads to various complications, therefore we briefly point out the method, the 
complications and how we have resolved them. 

The procedure starts with some initial configuration of link variables {U}\. We then 
update all links in lexicographic order, a process called a sweep, and arrive at a new 
configuration {t/} 2 . The updating process for each link proceeds as follows. Consider 
the link U{j. We identify which plaquettes contain this link: in three dimensions, there 
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are four such plaquettes. Now we calculate, for each element g G H, what the sum of 
the plaquette actions for each of these four plaquettes would be if were to have the 
value g. This gives a set of numbers 

where S 9k is the sum of the four plaquette actions with equal to g^. We now calculate 
a localized partition sum Z\j i : 

which can be used to calculate a set of probabilities {p{g)} g ^H for each group element g 



p{g) 



e Sg 



Z v ~ 

After a given number of sweeps no, the Monte Carlo algorithm arrives at the minimum 
of the action and the path integral expectation value of the operator O 

(o) = 1 jj^m (40) 

is given by taking the average of 0[{C/} n ], the value of O at gauge field configuration 

{U}n. 

^ no+m 

(0)MC estimate = ~ £ °^ U U (41) 

n=no+l 

However, for our purposes this scheme is troublesome for two reasons: it is tacitly 
assumed that the presence of the operator O in (40) does not change the value of the 
minimum of the action S and furthermore the loops of magnetic flux are very non-local 
objects and therefore highly unlikely to appear when using a local updating algorithm. 
This is illustrated in Figure [6] The shift upward of the functional S[U] is due to the 
presence of a magnetic flux string and the shift to the left is due to the non-locality 
of the magnetic excitations. The latter shift also occurs when a single loop of flux is 
inserted. 



The minimum of the action in the calculation of an S- matrix element (34) is altered 
by the insertions of the loop operators: the configuration for two non-commuting fluxes 
carries a string (see the discussion around Figure [5J that is massive and thus costs a 
finite amount of action. There is no way to get rid of this string and therefore the 
minimum value of the action in the presence of the two loops is shifted. We therefore 
have to amend the standard MC algorithm. Defining 

S = S m i n + 5S without operator insertion, 
S = + 5S with operator insertion, 

and noticing that around the minimum the actions behave identically, implying that SS 
and 8S are the same functions, expression (40) becomes 

f DUe-V™- s ^O[U]e- ss M 
\°f ~ fDUe- s m ■ (42) 
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S[{U}] 




{U} 

Figure 6. Schematic drawing of the action as a functional of the gauge field configuration {£/}. 
The insertion of non-commuting fluxes shifts the minimum of the action to a different location in 
the configuration space (due to the non-local nature of the excitations) and to a different value 
(due to the presence of a string). 

This leads to a modified Monte Carlo average 

. no+m 

(0)MC estimate = -e"^"^ £ 0[{U} n }. (43) 

n=no+l 

We now describe two approaches to the second problem in our MC measurements: 
the low probability that the local updating algorithm will converge to a gauge field 
configuration containing a (set of) magnetic flux loop(s). We will assume a single loop 
of pure magnetic flux is inserted, as nothing substantial will change in the case of 
multiple loops or the addition of dyonic charge. 

The first approach is based on the observation, illustrated in Figure [4], that we know 
the gauge field configuration (up to gauge transformations) that extremizes the action 
in the trivial vacuum with the insertion of a loop of magnetic flux: the /i-forest. We can 
therefore use this configuration as an ansatz in the MC algorithm. We start with a "cold 
lattice" , all links Uij = e, except for the /i-forest, for these links we set = h. This is 
an extremum of the action for the action if we set all @A^e — and /3 e ^> 1. To perform 
a measurement at some other value of the coupling constants, we can slowly change the 
coupling constants towards the desired values, performing a few MC updates after each 
step. The second approach is a more physical one. We initialize the lattice directly at 
the desired point in coupling constant space. The trick is then not to insert the loop 
all at once, but to slowly grow it, as illustrated in Figure [7j We start by twisting the 
action for four plaquettes around one link, as shown by the shaded plaquettes in the 
top left of Figure [7| After this, a number of MC updates are performed. Then the set 
of plaquettes that have a twisted action is changed as in the top right corner of the 
Figure. Again a number of MC updates is performed and so on. We have checked that 
in the trivial vacuum one obtains the /i-forest configuration using this procedure. Both 
of methods to insert flux loops have been used by us and we have verified that they lead 
to completely equivalent results. 
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Figure 7. Growing a flux loop in multiple steps. The shaded plaquettes have a twisted action, 
and the fat links show the convergence towards a /i-forest state. 

3.3. Results 

In this subsection we present the results of our Monte Carlo simulations. The first 
quantity we measured was the free energy as a means to map out a suitable subspace 
of the parameter space. It gives us an indication of the validity of our naive intuition 
about where nontrivial condensates should occur. 

Once we have found some region where symmetry breaking occurs we measure 
the open string expectation values to determine the respective condensates. After that 
we measure the unbroken and broken S"-matrix elements. Using the straightforward 
algorithm involving the auxiliary symmetries of our loop operators discussed in section 
2.3| allows us to find the branching matrix n® 1 as well as the S'-matrix of the effective 
U theory in the broken phase. 

Mapping out the phase diagram. The space of coupling constants in our theory is five- 
dimensional but it is not our goal to analyse it completely. We have restricted our 
search to some representative regions where nontrivial condensates do indeed occur. To 
study the location of the corresponding phase transitions we measured the free energy 



The modular S -matrix as order parameter for topological phase transitions 34 

F, which we define as the expectation value of the plaquette action, averaged over the 
spacetime lattice. The left plot of Figure M shows F as a function of (/3 e and (3 S ) and 












Figure 8. Plots of the free energy F for two-dimensional planes through the origin of the 
parameter space of the lattice model. In the left figure we have the (/3 e , /3 g ) plane and in the right 
figure we have the (/3 e = 0Xi > Pe) plane. See text for further comments. 

all other couplings equal zero. For small values of all the couplings appearing in the 



action (39), we are in the completely confining phase of the gauge theory, where all the 
open string operators of magnetic flux have a non-zero expectation value, and all loop 
operators carrying electric charge are confined. This corresponds to the plateau in the 
graph where F is maximal and tends to zero. 

The regions where the magnetic flux (e, 1) and (Xi,r°) have condensed can 
be anticipated on theoretical grounds by realizing that the coupling Pa is inversely 
proportional to the mass of flux A. In fact, when we look at the subgroup Ka generated 
by the elements in conjugacy class A, in particular 

K- e = {1,-1} 

K Xl = {1,-1,201,-20-!}, 

and set the couplings for the conjugacy classes containing the elements in Ka equal to 
one another, there is an extra gauge invariance U p — >• kU p for an element k G A in the 



plaquette action (|39j). In particular 

S = P(6 e (U p )+6e(U p )) 

is invariant with respect to U p — >■ — 1U P and 

S = p(8 e (U p )+8e(U p ) + 8 Xl (U p )) 

is invariant with respect to U p — >■ k U p , where k G {— l,ici, — i&i}- 

These left multiplications are exactly the kind appearing in the definition of the 



(loop) order parameters (33). Therefore one can establish, even without reverting to 
MC measurements, that the above actions, for large values of 0, produce the desired 
flux condensates. 
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(a) (Se(U p )) 



(b) (6e(U p )) 
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(d) <5x 2 (t/ p )) 



(e) (<M^)> 
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condensed 
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- U P =e Pe 

(f) Schematic of diagram 



Figure 9. Space-time avaraged expectation value of 5a{U p ) for each conjugacy class A. Shown 
is a (P e , /?e)-plane in coupling constant space where the other three couplings are zero. The color 
coding is such that red is the highest and blue the lowest value in each figure. In (f) we have 
identified the meaning of the various regions and the transition lines, where the red arrows indicate 
the trajectories used to determine whether the transitions are first or second order (see figure 12 
andpf)- 



One may verify this reasoning in the Figures [9] and 10 where we have probed the 
phase diagram more in detail by measuring the spacetime averaged expectation value 
of Sa(U p ) for all conjugacy classes A as a function of the relevant coupling parameters 
p. The red color indicates high values for the expectation value and we see that for 
all coupling parameters near zero all fluxes are condensed and thus all charges will be 
confined. This is what traditionally is called the "strong coupling phase (g ~ 3> 1)- 
Looking at the colorings for the various operators one readily identifies the various 
phases as indicated in the schematics of the subfigures (f). For example the symmetry 
with respect to the diagonal of the Figures [9^, and[9]D, shows that there are "Ising" like 
ordered phases, one with all plaquette values U p — e and the other with all U p = e. 
The in-between region is the region with the e flux condensate. Note that if the e flux 
would be the only one that phase would continue all the way to the origin, end we would 
exactly end up with the Z 2 pure gauge theopry phse diagram. 

In the region with e larger than approximately 2.0 and all other couplings near 
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(a) (S e (U p )) 



(b) (5e(U p )) 





(d) (S X2 (U P )) 



(e) (Sx 3 (U P )) 
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X i and e fluxes 

condensed 
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All fluxes 
condensed 

All charges 
confined 



(f) Schematic of diagram 



Figure 10. Space-time avaraged expectation value of 8a(U p ) for each conjugacy class A. Shown 
is a (/3 e , f3e — /?Xi)-plane in coupling constant space where the other two couplings are zero. The 
color coding is such that red is the highest and blue the lowest value in each figure. In (f) we 
have identified the meaning of the various regions and the transition lines, where the red arrow 
indicates the trajectory used to determine whether the transition is first or second order (see 



figure 12 : ). 



zero, the trivial vacuum is realized. All string operators with nontrivial magnetic flux 
have zero expectation value there. 



As pointed out in previous sections for example in relation (12), there is a very 



direct way to determine the condensate as well as the quantum embedding index q. 
This is by measuring the expectation value of the open string for each pure flux A and 
then summing over all fluxes. In the (e, 1) vacuum we obtain 



(M) 



(e,l) 



1.0, 



so q = 2, whereas in the (X\, T°) vacuum 



(e,l) 



1.0, 



l (*i,r°) 



2.0, 



so in this case q = 4. In Figure 11 we show the measurement of the vacuum expectation 
value for the (e, 1) open string as a function of the coupling constant (3 E , which 
demonstrates that such measurements clearly indicate where the transition takes place. 
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There is one more issue we like to address in our simulations, that is to determine 
the order of the transitions we have identified. A conventional approach is to search 
for a hysteresis effect across a first order transition, but because of the relative modest 
size of the lattices used this is not an optimal approach. A method that is working 
much better is to directly probe the system at a given sequence of coupling constants 
around the transition and to see whether there is a coexistence region where both phases 
occur in the samplingFj To perform these measurements we use the parallel tempering 
method [39] to overcome local minima in the action landscape. The idea behind this 
method is to initialize a range of lattices simultaneously, all at different couplings along 
a trajectory in coupling constant space starting in phase one and ending in phase two. 
The updates of this ensemble then consist of the updates of each of the individual lattices 
and, occasionaly, a swap of two adjacent lattices. The swap between lattices 1 and 2 is 
accepted with a probability 

exp(Si(l) + S 2 (2)) 



p(l <+2) =min^ 1, 

I exp (S 2 (l) + <Si(2)) 

where Si (2) means using the action (in particular, the set of couplings) of lattice 1 to 
evaluate the field configuration of lattice 2 et cetera. One can prove that this satisfies 
detailed balance. In effect, each lattice will perform a random walk through coupling 
constant space along the chosen trajectory, allowing a "cold" lattice to thermalize in 
the "high temperature" region, thus overcoming the local minima of the action. 

We have made measurements for the trajectories indicated by the arrows in the 



Figures p5f and 10 :. The results of these measurements for the horizontal arrow is given 



in Figure 12 and for the vertical arrow in Figure 13 We find that in that the horizontal 
trajectory the transition from the strongly coupled phase corresponding to the left peak 
in Figure 12 to the trivial phase corresponding to the right peak indeed goes through 
a coexistence region corresponding to the values of the coupling parameter where both 
peaks are present as in subfigures [T2|d and [life. 

The result for the vertical trajectory corresponding to the transition from the trivial 
phase to the broken (Xi,r ) phase is given in figure 13, where we see that the peak 
shifts continuously implying that the transition is second order. We can understand 
this transition as follows. In this region of coupling constant space, all fluxes except 
the e flux are very heavy. This means the ground state is essentially that of a Z 2 
gauge theory. Since Z 2 gauge theory in 3 dimensions is Kramers- Wannier dual to the 3 
dimensional Ising model, it has the same phase structure [31] . We therefore expect this 
phase transition to lie in the same universality class. 



Measuring the (broken) modular S-matrices. We have measured the (broken) S-matrix 
elements using the simple algorithm involving the auxiliary symmetries of our loop 
operators. This allows us to obtain the unbroken S-matrix as well as the branching 

+ We would like to thank Simon Trebst for pointing this out to us. 
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Figure 11. Vacuum expectation value of the (e, 1) open string as a function of the coupling 
constant /3 g , showing that the nonlocal open string operators are good order parameters to 
characterize topological phase transitions. Length of the string: 4 plaquettes, measurements 
on a 16 3 lattice, (3 e — 3.0, other couplings zero. 



matrix n" 1 and the S-matrix of the effective IA theory in the various broken phases. 
Here we exploit the relation (27) for the measurement, and relation (18) : 

a it bj 

relating S uv to the measured S'-matrix in the broken phase. We first measured the 
unbroken S'-matrix in the D(D 2 ) phase and obtain results identical to the matrix 
calculated using defining formula (|6| , the result is given in Table [2] and is of course 
also consistent to the matrix obtained from the relation (18) with $ = 0. The accuracy 
of the measured matrix elements in represented in the table as integers is smaller than 
5%. 

The branching matrices n a u can be obtained from measuring the broken S-matrices. 
The columns in these matrices correspond to the different U sectors. If we see two rows 
or columns with different parents a, b in the A theory that are proportional to each 
other, a and b branch to the same U sector u. Conversely, if different u fields correspond 
to the same a field that means that the a splits in the broken phase. We have listed the 
results for the broken S-matrix in the (e, 1) vacuum in Table [3] 

To realize the splittings between the irreducible representations using the auxiliary 
gauge symmetries alluded to in section 2J3 we found the following construction to suffice, 
(e, 1) vacuum 

• (X^T )! is realized by the operator A^ Xi,r °\ 



The modular S -matrix as order parameter for topological phase transitions 39 



(1.906, 0.0. 0.0, 0.0, 0.0) (1.916, 0.0, 0.0, 0.0, 0.0) (1.926, 0.0, 0.0, 0.0. 0.0) 

N( <,B f {U p )} ) N( (6,(U P )) ) N( (6 r {<J P )) ) 




0.5 0.6 0.7 0.8 0.9 1.0 0.5 0.6 0.7 0.8 0.9 1.0 0.5 0.6 0.7 0.8 0.9 1.0 



(a) (3 e = 1.906 (b) f3 e = 1.916 (c) /3 e = 1.926 



(1.936, 0.0, 0.0, 0.0, 0.0) 

N( <&,((/,)>) 




(d) f3 e = 1.936 

Figure 12. The sequence of plots is across the transition from the strongly coupled phase with all 
fluxes condensed and all charges coinfined, to the trivial phase. This trajectory corresponds to the 
horizontal arrow in figures [9J" and[lOf, where 1.906 < f3 e < 1.936 and all other couplings equal 
zero. Plotted along the x-axis is average expectation value of the percentage of trivial plaquettes 
with Up = e and along the y-axis we plot the number of times that that percentage is measured in 
a simulation of 4000 runs on a 10 3 lattice. The figures clearly show a shift from peak on the left 
to on the right, with a double peak in between, this is the signature of region where both phases 
coexist, i.e. of a first order transition. 
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Table 2. The S'-matrix for the (unbroken) D(D2) theory (up to the normalisation factor 
1/Da = 1/8) as measured in the trivial vacuum. We put integers in the table as the accuracy is 
below the 5%, i.e. 1 actually stands for read as 1. ± 0.05. 

• (Xj,r°) 2 is realized by the operators A( Xi,r °)A( e ' J ^. 

• (Xi,r 2 )i is realized by the operator A^ Xl ' r2 ^ with {x^ = e,x_j CT1 = icr 2 }. 
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Figure 13. The sequence of plots is across the transition from the trivial phase to the phase 
with the e condensate, the trajectory corresponds to the vertical arrow in figure [9f, where 
1.42 < fa < 1.72, fa — 3.0 and all other couplings equal zero. Plotted along the x-axis 
is average expectation value of the percentage of trivial plaquettes with U p = e and along the 
y-axis we plot the number of times that that percentage is measured in a simulation of 4000 runs 
on a 10 3 lattice. The figures only feature only a single peak that smoothly moves from one phase 
to the other, indicating a smooth second order transition, presumably corresponding to the 3-d 
ising model transition. 



• (X!,r 2 ) 2 is realized by the operators A^ Xl,r2 "> A^ 1 ^ with {x iai = e, x_ ioi = ier 2 }. 

• (X 2 ,r 2 )! is realized by the operator A^ 2 ' r2 ) with {x i(T2 = e,x_ ia2 = io~i}. 

• (X 2 ,r 2 ) 2 is realized by the operators A^ 2 ^ A^' 1 ^ with {xi U2 = e,X-i a2 = iai}. 

• (X 3 ,r 2 )i is realized by the operator A( X3 ' r2 ) with {xi a3 = e,X-i a3 = io~i}. 

• (X 3 ,r 2 ) 2 is realized by the operators A( X3 ' r2 )A( e ' j3 ) with {xi U3 = e,X-i a3 = iai}. 

(X\, r°) vacuum 

• (Xj,r°)i is realized by the operator A( Xi,r °). 

• (Xi,r°) 2 is realized by the operators A^ 1 ' r °)A^ Jl ). 

• (X 2 ,T 2 )i is realized by the operator A( X2 ' r2 ^ with {xi a2 = e,X-i a2 = iai}. 

• (X 3 ,r 2 )i is realized by the operator A^ 3 ' r2 ^ with {xi a3 = e,X-i a . 3 = io\}. 

We see that in Table [3] the columns (rows) for the sectors (e, a) and (e, a) for 
a = 1, Ji, J 2 , J 3 are identical and thus that the corresponding fields have to be identified. 
This leaves us with 16 sectors for the broken U theory. Summing the entries as prescribed 
by formula (3.3) yields exactly the S'-matrix of the /J(Z 2 £g> Z 2 ) theory, which is given 
in Table H 



In Table[5]we have listed the result for broken S'-matrix in the {Xi, T°) vacuum, here 
we have to identify the sectors (e, 1), (e, 1) and (X 1 , T°) 1 , the sectors (e, Ji), (e, J\) and 
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Table 3. The broken S-matrix as measured in the (e, 1) vacuum, where the cloumns and rows 
of zeroes corresponding to the confined fields are left out. Identifying identical columns and rows 
we obtain the familiar S-matrix of the D(Z 2 <S> Z 2 ) theory. 
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Table 4. The S'-matrix of the U — Di^L^®^) theory (up to normalisation factor of 1/Djj = 1/4) 
as obtained from the broken S'-matrix measured in the (e, 1) vacuum, after leaving out the rows and 
columns with only zeroes of the confined fields and after identifying identical rows and columns. 



(X!,r ) 2 , the sectors (X 2 ,r°)i and (X 3 ,T ) h and the sectors (X 2 ,r 2 )i and (X 3 ,r 2 )i. 



These results are all fully consistent with the algebraic analysis presented in Section 3.1 



Let us illustrate the method by calculating a few sample S'-matrix elements in the 
(Xi,r°) condensed vacuum. The U theory should be D(Z 2 ); let us first calculate the 
S'( +i+ )( +i+ ) element, the (+, +) sector being the new vacuum 



Q 



(%,l)(e,l)> $ + (%l)(e,l)} $ + (%l)(X 1) rO)> $ + 

(S(x lt ro)( e ,i))^ + (S( Xl ,ro)(e,i))^ + vS(x 1) ro)(x 1 ,ro)) 4 

-^(1 + 1 + 2 + 1 + 1 + 2 + 2 + 2 + 4) = l, 
48 2 
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Table 5. The broken S'-matrix for the D{D2) theory as measured in the (Xi,r°) vacuum. 

in agreement with Table [6j The contributions to the above matrix element would be 
equal if we had used the S'-matrix elements as measured in the trivial vacuum. 























o 
U 


u 












CI 






oT 




* 




0,i) 


(+,+) 


1 


1 


1 


i 


0, Ji) 


(+, -) 


1 


1 


-1 


-l 


(x 2 ,r°) 1 


(-,+) 


1 


-1 


1 


-i 


(* 2 ,r 2 )i 


(-,-) 


1 


-1 


-1 


i 



Table 6. The modular S'-matrix for the D(7L2) theory (up to normalisation factor 1/Djj = 1/2) 

To appreciate the importance of the measurements in the broken vacuum, consider 
the matrix element S , (_ + )(_ _). The parents of the (— ,+) sector are (X 2 ,r°)i and 
(X 3 ,T°) 1 and those of the (-, -) are (X 2 ,r 2 )! and (X 3 ,r 2 )!. 




~((-4) + (-4) + (-4) + (-4)) = -i 

. We see that after completing the calculation along this line we obtain the S'-matrix 
of the .D(Z) theory, as given in Table [6j Note that if we had used the S- matrix of the 
unbroken theory, the S(x 2 ,r°)(x 3 ,r 2 ) and S(x 3 ,r°){x 2 ,r 2 ) would have been zero. 

4. Conclusions and outlook 

In this article we have studied euclidean lattice models for Discrete Gauge Theories. 
We have introduced a set of multiparameter actions for these theories that display a 
rich phase structure, and showed in particular that all the allowed condensates of pure 
magnetic flux are realized in certain well anticipated regions of coupling constant space. 
The set of open string operators that we defined form a set of order parameters that 
allowed us to determine the content of the condensate and to measure the topological 
symmetry breaking index q. 
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Once the condensate is identified, we have shown how to unambiguously reconstruct 
the S-matrix of the low-energy theory in a broken or unbroken phase by measurements 
of the complete set of braided loop operators, using the anyonic loop operators we 
proposed in earlier work [16]. Due to an auxiliary gauge symmetry these operators are 
particularly well suited to detect the nontrivial splittings of fields that correspond to 
fixed points under fusion with the condensate. We found that as expected the excitations 
that are confined in a broken vacuum give rise to rows and columns of zeroes in the 
broken S'-matrix. Our work clearly demonstrates that the euclidean approach allows for 
a very straightforward method to completely determine the nature of the broken phase. 

Our work showed that the reason the modular S'-matrix changes in the broken 
phase is largely due to the contribution of the so-called vacuum exchange diagram. In 
an upcoming more theoretical paper [2FJ we will extend the approach used in this work, 
the use of observables and in particular the S-matrix to determine the phase structure 
of a TQFT to a far wider range of theories, in particular the SU(N)k TQFT arising 
from Chern-Simons actions. 

It would be interesting to study different models exhibiting different topological 
phases by somehow formulating them in the euclidean 3-dimensional framework, to our 
knowledge such an approach is unfortunately not yet available for Chern Simons theories. 
One expects that for Levin Wen models [ID] our approach could be implemented though. 
Another path is to investigate the phase structure after adding dynamical matter fields 
that transform nontrivially. It is known that in such situations the Wilson type criteria 
break down as the strings can break, this necessitates the development of different 
diagnostic tools [351 E] • 

The authors would like to thank Jan Smit, Joost Slingerland, Simon Trebst and 
Sebas Eliens for useful discussions. JCR is financially supported by a grant from FOM. 
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Appendix A. Fusion rules for D 2 DGT 
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